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ABSTRACT 

We present an overview of a 90-orbit Hubble Space Telescope treasury program to obtain near 
ultraviolet imaging of the Hubble Ultra Deep Field using the Wide Field Camera 3 UVIS detector 
with the F225W, F275W, and F336W filters. This survey is designed to: (i) Investigate the episode 
of peak star formation activity in galaxies at 1 < z < 2.5; (ii) Probe the evolution of massive galaxies 
by resolving sub-galactic units (clumps); (iii) Examine the escape fraction of ionizing radiation from 
galaxies at z ~ 2 — 3; (iv) Greatly improve the reliability of photometric redshift estimates; and 
(v) Measure the star formation rate efficiency of neutral atomic-dominated hydrogen gas at z ~ 1 — 3. 
In this overview paper, we describe the survey details and data reduction challenges, including both 
the necessity of specialized calibrations and the effects of charge transfer inefficiency. We provide a 
stark demonstration of the effects of charge transfer inefficiency on resultant data products, which 
when uncorrected, result in uncertain photometry, elongation of morphology in the readout direction, 
and loss of faint sources far from the readout. We agree with the STScI recommendation that future 
UVIS observations that require very sensitive measurements use the instrument's capability to add 
background light through a "post-flash" . Preliminary results on number counts of UV-selected galaxies 
and morphology of galaxies at z~l are presented. We find that the number density of UV dropouts 
at redshifts 1.7, 2.1, and 2.7 is largely consistent with the number predicted by published luminosity 
functions. We also confirm that the image mosaics have sufficient sensitivity and resolution to support 
the analysis of the evolution of star-forming clumps, reaching 28-29th magnitude depth at 5<r in a 0'.'2 
radius aperture depending on filter and observing epoch. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: high-redshift — 
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1. INTRODUCTION 

The great succ ess of the GALEX mission 
(Thilkcr ct al.l 120051 ) revolutionized the study of 
galaxies in the ultraviolet (UV). But it has left us in 
the curious position of having extraordinary detail on 
the UV emission and structure of the closest galaxies 
(from GALEX) and quite distant ones (where the UV 
rcdshifts into optical bands), but having significantly 
less data in between. The rest-frame 1500 A continuum 
(FUV) is an important tracer of star-formation, because 
it samples the output from hot stars directly. The 
star-formation density of the Universe peaks in the 
epoch 1 < z < 3, which requires deep near-ultraviolet 
(NUV; A ~ 2000 - 3500 A) observations to measure the 
redshifted FUV. 

A new generation of Hubble Space Telescope (HST) 
surveys have been approved to begin filling this gap 
through deep, high spatial resolution imaging. The Wide 
Field Camera 3 (WFC3) UVIS channel provides revo- 
lutionary sensitivity in the NUV. Shortly after instal- 
lation, the WFC3 team conducted Early Rel ease Sci- 
ence observations fERS: IWindhorst et al.ll20Tll ). includ- 
ing a first look, multi-wavelength extragalactic survey. 
The ERS included about 50 square arcminutes of NUV 
imaging, at 2,2,1 orbit depths in the F225W, F275W, 
and F336W filters respectively, reaching 26.9 magnitudes 
(AB). More recently, the Cosmic Assembly Near-IR Deep 
Extragalactic Legacy Survey (CANDELS; iGrogin et al.1 
120111 : iKoekemoer et al.l 1201 lh has started observations 
with UVIS. CANDELS will observe the northern field of 
the Great Observator ies Origins Deep Survey(GOODS 
iGiavalisco et al.ll2004D with the F275W filter in the con- 
tinuous viewing zone, for a total predicted depth of 27.0 
magnitudes (AB; 5a in a 0'.'2 radius aperture) over 78 
square arcminutes. 

In this paper, we describe a new program (GO- 
12534; PI=Teplitz) to obtain deep, N UV imaging of 
the Hubble Ultradeep Field (UDF; iBeckwith etaTl 
2006). The UDF provides one of the most stud- 
ied fields with a wealth of multi- wavelength data 
dRosati et al.l 120021: IPfizkal et al.l[200l [Yan et all 120041: 



Thompson et al. 2005; Labbc ct al. 2006; Kaiisawa ct al 
2006 1 : IBouwens et al.ll2006l: lOesch et al.ll2007tlSiana et alj 
20071: iRafelski et all [20091: 
Vover et al.1 120091: iRe tzlaff et al 



Nonino et al.l 120091: 
20101: IGrogin et all 



20111: IKoekemoer et all 120111 : IBouwens et all 12011 



Elbaz et all 120111: iTeplitz et all 120111 : IKoekemoer et all 
2012tlEllis et al.M2013h . enabling the best return on this 
new investment of telescope time. This project obtained 
deep images of the UDF in the F225W, F275W, and 
F336W filters at 30 orbit depth per filter (see Figure 
[J), with the goal of reaching 28-29th magnitude (AB) 
depth at 5a in a Of.' 2 radius aperture. The program 
was designed to use 2x2 onboard binning of the CCD 
readout to improve sensitivity. That mode was only 
used for the first half of the observations, at which point 
it became clear that another strategy is better. The 
second half of the observations were obtained without 
binning of the CCD readout, but with the UVIS ca- 
pability to add internal background light, "post-flash", 
to mitigate the effects of degradation of the charge 
transfer efficiency of the detectors. We will discuss the 
implications of these choice for both sensitivity and 
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Fig. 1.— Throughput of the WFC3-UVIS filters used in the 
UVUDF: F225W in blue, F275W in green, and F336W in red. 
These throughputs include the quantum efficiency of the CCD. 

data reduction. Combined wit h previous imaging of the 
UDF in the far- ultraviolet (see lSiana et alll2007|) . these 
new observations (hereafter UVUDF) will complete the 
pan-chromatic legacy of this deep field. 

We describe the science goals of the project in Section 
[21 survey strategy and observations in Section [3] we out- 
line data reduction and source extraction in Section^ we 
characterize the data quality and discuss issues related 
to the charge transfer efficiency of the CCD in Section 
El In Section [6] we describe preliminary analysis of the 
data and initial science results, before summarizing in 
Section [71 Throughout, we assume a A-dominated flat 
universe, with Ho = 71 km s _1 Mpc -1 , 51a = 0.73, and 
n m = 0.27. 

2. SCIENCE GOALS 
2.1. Tracing the evolution of star formation 

Observations of UV-bright galaxies trace the evolution 
of cosmic star formation and provide key constraints on 
galaxy formation. The UVUDF detects galaxies with 
star formation rates (SFRs) greater than ~ 0.05 M Q /yr 
at z ~ 2 — 3 (in the absence of dust extinction) with the 
same UV color selection techniques used at higher red- 
shift. For example, the Lyman break galaxy (LBG) selec- 
tion, whereby high redshift galaxies are identified by their 
strong flux decrement at short wavelengths due to the Ly- 
man break fe ature, is routinely used in man y studies (e.g. 
Steidel et all 119991: lAdelberger et all 120041: iReddv et all 
20081 : IBouwens et al.ll2011l) . When more photometric in- 
formation is available, more complex methods become 
available (see Section 12. 4|) . Measuring the combination 
of the UV luminosity function and the mass function of 
UV-selected galaxies will provide a statistical picture of 
the history of star formation in these sources, i n red- 
shift slices between 1 < z < 2.5 (|Lee et al.ll2012bO . UV- 
selection in this epoch will enable significant spectro- 
scopic follow-up, with access to vital rest-frame optical 
diagnostics of extinction, metallicity, and feedback. We 
provide an initial LBG selection in UVUDF in Sectional 

One of the largest sources of systematic error in esti- 
mates of the SFR and the cosmic star formation history is 
the fact that dust absorbs and reproce sses approximatel y 
half of the starlight in the universe ((Kennicuta [l998a) . 
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The amount of re-radiated light, quantified by the ra- 
tio of integrated IR to UV luminosity, IRX=L/fl/L[/y, 
has been found to be correlated with the UV spectral 
inde x, (where fx oc A ^), in local starburst galaxies 
(e.g. IMeurer et all Il999l ). This correlation is routinely 
used to correct UV SFR estimates for dust attenua- 
tion in highly st ar forming galaxies at all redshifts (e.g. 
LBGs and BzKs:IAdelberger fe Steidell2000HDaddi et all 
20071: iReddv et all 120101 12012at iKurczvnski et al.ll2012t 
Lee et afll2012al ). UV bright galaxies and IR luminous 
galaxies (Lib, > 10 11 L Q ) at lower redshifts are found to 
be br oadly consistent wit h the starburst IRX-/3 correla- 
tion (jOverzier et al.ll20TTI ). Understanding the effects of 
extinction at high redshift requires detailed study of nor- 
mal galaxies 7-10 Gyr into the past (the epoch probed 
by the UVUDF), where both the UV slo pe and the in- 
frared emission can be measured, (e.g. iBoissier et al l 
20071 iSiana et al.ll2009l: [Swinbank et al.ll2009t iBuat et all 



2010; Bo uwens et al.ll2012tTFinkelstein et al.M 2012). 



2.2. The Build-up of Normal Galaxies 

The role (and nature) of feedback, and the relative im- 
portance of merging in galaxy mass growth are still de- 
bated issues. Observations show that "normal" galaxies 
were in place at z ~ 1, with stellar population and scal- 
ing relations consistent with passive evolution into the 
hom ogeneous population observed in the loca l Universe 
(e.g. iScarlata et all 120071: iCimatti et alll2008t ). This sit- 
uation changes drastically looking back just a few Gyrs. 
Among the diversity and complexity of massive galaxy 
types, two types have been extensively studied: gas- 
rich clumpy disks forming stars at rates of 100 M /yr 
that do not have counterparts in the local Universe (e.g. 
Daddi et all 120101: lElmegreen et all 120051: iGenzel et all 
20081 ). and passive objects that are obser ved to be ~30 



times denser than any galaxy today (e.g. ICimatti et all 
120081: Ivan Dokkum et alJl2008h . The former are key to 
understand the role of instability and gas accretion in the 
formation of disks and bulges (by migration and merging 
of the clumps) ; the latter because we do not yet under- 
stand the physics of quenching of star formation and the 
role that compactness plays in it. 

It is tempting to think of these well-studied popula- 
tions as different phases in the formation of local galaxies. 
Secular evolution of star-forming sub-structures within 
gas-rich disks could lead to the formation of bulges, 
and the compactness of high— z spheroids would be the 
result of the highly d i ssipative merger of the clumps 
(jElmegreen et all 120081 : iDekel et all 12009ft . Clumps are 
predicted to be fueled by cold (T < 10 4 K) gas streams 
that efficient ly penetrate the hot medium of the dark 
matter halo (jKeres et al.ll2009| ). The UV morphology of 
LBGs at z = 3 — 4 are a lso suggestive of this process 
([Ravindranath et alll2006T) . Furthermore, it is still not 
clear what mechanism quenches the star formation in the 
newly formed bulges, what prevents more gas from cool- 
ing and forming stars, and what drives the size evolution 
of compact spheroids. 

Current HST observations allow us to derive stellar 
masses, SFR, surface density of star formation, and the 
extinction of individual bright clumps at z ~ 2 — 3 by 
fitting the spectral energy distribution (SED). However, 
without access to the rest-frame UV, our assessment 



of star-formation activity becomes poorer at lower red- 
shifts. At z ~ 2.3, such structures are found to have sizes 
of ~ 1.8 kpc, typical masses of several 10 7 M fn, and ages 
of ~ 0.3 Gyr (jElmegreen fe Elmegreenll2005ft . 

The UVUDF observations are designed to provide the 
depth and resolution (~ 700 pc) to study sub-galactic 
structures at 0.5 < z < 1.5 at consistent rest-frame UV 
wavelengths, offering continuity with measurements at 
low and high redshift. We confirm the utility of the data 
for this purpose in Section [6] Measurement of the typ- 
ical UV sizes and luminosity will constrain stellar-mass 
and stellar-population properties using the full SED. Fi- 
nally, the data will enable comparison of the colors of 
individual sub-galactic units at different radii for the SF 
galaxies at z < 2 and z — 3. A color gradient would 
be expected if there is migration of previously formed 
structures towards the center to form the bulge. 

2.3. Contribution of galaxies to the ionizing background 
(below 912 A) 

Star-forming galaxies are likely responsible for reion- 
izing the universe by z ~ 6, assuming that a high frac- 
tion of Hi-ionizing (Lyman continuum; LyC) photons 
are able to escape into the IGM. Recent studies sug- 
gest t hat the escape fract i on, f eKr is higher at high red- 
shift dShaplev et al.ll2006t llwata et al . 2009; Sian a et "all 



20101: iNestor et all 120131 but see 



2010t Bridge et al 

Vanzclla et al. 02011 . However, it is extremely difficult 
to directly measure the LyC at z > 4 due to the increas- 
ing opacity of the IGM. Thus, it is important to under- 
stand the physical conditions that allow LyC escape at 
2 < z < 3 and to determine if those conditions are more 
prevalent during the epoch of reionization. 

Ground-based surveys suffer from significant fore- 
ground contamination, and from not knowing from which 
part of the source the ionizing emission is escaping. HST 
resolved images of both the ionizing and non-ionizing 
emission of galaxies are necessary to confirm the ex- 
treme ionizing emissivities suggested b y previous surveys 
(|Iwata et all 120091 INestor et al.l[2013h . The UVUDF fil- 
ters will enable measurement of the LyC escape frac- 
tion at redshifts z ~ 2.20,2.45,3.1 in F225W, F275W, 
F336W, respectively (see Figure Q] for the filter through- 
puts). 

2.4. Improved Photometric Redshifts 
Despite intensive spectroscopic surveys that have 



provided hund re ds o f redshifts dSzokoly et all 



2005 



2004; 



2006 



2010; 



| Le Feyre et all 120051: iVanzella et al . 
20081 I2009t iPopesso et all 120091: iBalestra et all 
iKurk et al.l 120131 ). the majority of sources in the UDF 
are either too faint or otherwise inaccessible. Redshifts 
must therefore be determined either through color 
selection or photometric redshifts (photo-z). However, 
young star-forming galaxies often lack strong continuum 
breaks in the rest-frame optical, making accurate 
photo-zs nearly impossible with only optical+NIR data. 

The three UVIS filters target the dominant signature 
of the galaxies' SEDs in the redshift range 1.2 < z < 
2.7 - the Lyman break. This feature will allow color 
selection of these galaxies, and will resolve many of the 
photo-z degeneracies and thereby improve the photo-z 
fits. While photo-z's currently exist for all objects in the 
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UDF (|Coe et al.1 120061) . they often have multiple peaks 
in their probability distribution fu nctions, P(z), making 
the true redshift uncertain. In fact, Rafclski et al. (2009) 
found that the introduction of the ground-based u-band 
data improved the photo-z's for 50% of the z ~ 3 sample. 
However, their results suffered from the limited angular 
reso lution and depth of ground-based u-band data (see 
also lNonino et all 12009). The F336W filter significantly 
improves the redshifts from 2 < z < 3 and z < 0.2, while 
the F275W filter improves the redshifts at 1.5 < z < 2 
and z < 0.3, and the F225W filter improves them at 
1 ^ z < 1.5. 

2.5. Star Formation Rate Efficiency of Neutral 
Atomic- Dominated Hydrogen Gas 

Th e locally established Kennicutt-Sc hmidt (KS) rela- 
tion (jKennicuttl I1998U iSchmidtl I1959D relates the gas 
density and the SFR per unit area, Hsfr °c ^gas- While 
this assumption is reasonable at low redshift for typi- 
cal galaxies, it has been shown not to h old for neutral 
atom i c-dominated hydrogen gas at z ~ 3 ( W olfe &: Chen! 
120061: iRafelski et ail 120111 ). Nonetheless, cosmologi- 
cal simulations often use the KS relation at all red- 
shifts, for both at omic and molecular hydrogen gas (e.g. 
iKeres et al.l[2009l ). 

Damped Lya systems (DLAs; see IWolfe et al.l[2005l for 
a review) are large reservoirs of neutral hydrogen gas. 
At z ~ 3, the in situ SFR of DLAs is found to be 
less than 5% of wha t is expected from the KS relation 
(|Wolfe fc Ched 120061) . This means that a lower level of 
star formation occurs in DLAs at z ~ 3 than in modern 
galaxies. Another possibility is that in situ star forma- 
tion may occur at the KS rate only in DLA gas associated 
specifically with LBGs. These DLAs are constrained by 
measuring the spatially extend ed low-surface - bright ness 
(LSB) emission around LBGs. IRafelski et aD (|2011l ) de- 
tect such emission on scales up to ~ 10 k pc in a sam- 
ple of z ~ 3 LBG's (|Rafelski et alJl20M ) in the UDF 
F606W image (rest-frame FUV). The emission is mea- 
sured to > 31 mag arcsec -2 and on large scales by stack- 
ing z ~ 3 LBGs that are isolated, compact, and sym- 
metric. The resulting SFR around LBGs was found to 
be ~ 2 — 10% of what is ex pected from the local KS 
relation (jRafelski et al.ll2011l) . 

This result can b e used to constrain models of galaxy 
formation at z ~ 3. iGnedin &: Kravtsovl <|2010f ) conclude 
that the main reason for the decreased efficiency of star 
formation is that the diffuse ISM in high redshift galaxies 
contains less dust, and therefore have a lower metallicity 
and a lower dust-to-gas ratio, which is needed for shield- 
ing in order to cool the gas and form stars. This notion 
matches the o bservation that DLA metallicities decrease 
with redshift {Rafclski et al.ll2012l ). and therefore we ex- 
pect that the efficiency of star formation may be corre- 
lated with redshift. This effect must be further under- 
stood and taken into account when interpreting models 
of galaxy formation and evolution. 

The transition from the lower star formation efficien- 
cies at z ~ 3 to those on the Hubble sequence at z ~ 
may be observable at redshifts in between. We plan to 
find that transition or constrain when and how it oc- 
curs by probing the star formation in the LSB regions 
around moderate redshift LBGs. It is only in the outer 



diffuse regions, where the metallicity is lower, that the 
KS relation is seen to be evolving. The NUV coverage 
of the UDF enables us to detect this star formation at a 
range of intervening redshifts by providing significantly 
improved photo-z's (Section 12. 4p at z ~ 2 — 3 in order 
to identify LBGs to stack in the optical UDF data, and 
possibly by stacking the UV data themselves at z ~ 1, if 
the CTE corrected data permit (see Section 15.1.11) . 

3. OBSERVATIONS 

The UVUDF program was executed in three epochs, 
due to the heavy scheduling constraints on HST in Cy- 
cle 19 (Fall of 2011 through Fall of 2012). Table Q] lists 
the orbit distribution and position angle of each set of 
observations. In each case, a common pointing center 
is used: RA: 03 32 38.5471 DEC: -27 46 59.00 (J2000). 
Figure [2] shows the orientation of each epoch compared 
to previous UDF programs. 

The UVIS focal plane consists of two CCDs, each with 
4146 x 2051 pixels. The plate scale is 0'.'0396/pixel at 
the central reference pixel. After accounting for the over- 
scan regions, the usable area of each CCD is 4096 x 2051 
pixels. There is a physical gap between the CCDs that 
corresponds to about 30 pixels (1'.'2). 

WFC3/UVIS has a field of view of 162" x 162", larger 
than the WFC3/IR channel (136" x 123") but smaller 
than the optical field of the Advanced Camera for Sur- 
veys' Wide Fiel d Camera (ACS/WFC; 202" x 202"; 
iFord et alJ l2002T >. The UVUDF observations are well 
matched to the WFC3/IR pointings from two observ- 
ing programs, as shown in Figure [3J The first pro- 
gram (GO-11563, PI=Illingworth) was excecuted in 2009 




(HUDF09 & 12) , 

Fig. 2. — The footprint of the the UVIS pointing for epochs 1, 
2, and 3 are shown as blue squares, with each epoch individually 
labeled. The greyscale image is the V-band ACS image of the 
UDF from Bec kwith eTaLl (2009 ). with North up and East left. 
The shaded regions are the footprints of ot her HST imaging of th e 
UDF. The orange represent s NICMOS IR IIThompson et al .([20051 1 . 
the green ACS-SBC FUV l|Siana et a l. 2007), an d the red WFC 
near-infrared imag ing from HUDF09 and HUDF12 (Bouwc ns et al.l 
2011; Ellis et al. 2013). The readout direction is perpendicular and 
away from the blue lines marking the chip gap in each epoch, such 
that the readout is located furthest from the chip gap. 
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TABLE 1 
UVUDF Observing Epochs 



Epoch 


Observing Window 


ORIENT 1 


Orbits per UVIS filter 


Orbits per ACS filter 


Epoch 1 


March 2 -March 11 2012 


96.0 


6 


4:3:ll a 


Epoch 2 


May 28- June 04 2012 


197.25 


8 b 


20:2:2:2 C 


Epoch 3 


August 3-September 7 2012 


264.0 


16 d 


46 



Note. — List of orbit distribution and position angle for each set of observations. 
1 The ORIENT keyword is denned in Section [Q 
a Parallel orbits per filter in order F435W:F606W:F814W. 

b Due to two failed visits, F336W has 10 orbits per filter, while F275W and F225W have 8. 
c Parallel orbits per filter in order F435W:F606W:F775W:F850L. 
d F336W has 14 orbits per filter, while F275W and F225W have 16. 



(HUDF09; IQesch et all [2010cllbl : IBouwens et all I20T1 . 

The second program (GO-12498, PI=Ellis) was executed 
after UVUDF a t the same pointing as the HUDF09 
(|Ellis et al.N2013f) . The footprint of previous UV imag- 
ing of the UDF taken with the ACS Solar Blind Chan- 
nel (SBC) (JSiana et al.H2007r> and IR imaging taken with 
NICMOS llThoiiipMm et all I2005D is also shown in the 
figure [2] 

Observations were obtained in visits of 2 orbit dura- 
tion in order to maximize schedulability. Each visit used 
a single UVIS filter. These visits were linked in groups 
of 3 in the scheduling instructions to guarantee that all 
three filters were obtained at the same oriention. During 
each 2-orbit visit, four exposures were taken. Typically 
this schedule allowed about 1300 seconds of integration 
per exposure. In total, we obtained ~ 82,000 seconds of 
integration per filter in the full overlap region (see Ta- 
ble [2]). Half the data were taken with binning of the 
CCD readout, while the other half were taken without 
binning, but with the use of the post-flash capability 
(see Section [3. lj) . The unbinned Epoch 3 exposures were 
dithered with the standard WFC3-UVIS-DITHER-BOX, 
which is a 4 point dither pattern with a point spacing of 
0'.'173. The binned Epoch 1 and 2 exposures are dithered 
in an analogous way, but with doubled spacing of 0'.'346. 

An exception to the observing plan occurred in two 
visits ("IN" and "2H" in the HST schedule), resulting in 
the loss of both visits in Epoch 2, one for F275W and one 
for F225W. These visits were rescheduled during Epoch 3 
(as visits "5N" and "6H"), and executed as planned at 
that time. 

The area of full overlap between dithered exposures, 



and thus full sensitivity, is 6.2 arcmin , or 



of the 



area of the UVIS detector. The full NUV UVIS overlap 
region and all of Epoch 3 are completely covered by the 
deep ACS optical data. The footprint of the UVIS point- 
ing is overlaid on the ACS F606W image of the UDF in 
Figure H The full WFC3/IR pointings (HUDF09 and 
HUDF12) are covered by the NUV UVIS data. 

3.1. Charge Transfer Inefficiency 

Over time, radiation causes permanent damage to the 
CCD lattice, decreasing the charge transfer efficiency 
(CTE) during readout. The CTE degradation is a seri- 
ous problem for low background imaging of faint sources, 
resulting in decreased sensitivity and uncertain calibra- 
tion for extended sources. The effect is worse for ob- 
jects that are far from the CCD readout, that is for ob- 
jects close to the gap between the two detectors in the 
case of UVIS. The degradation of the UVIS CCDs has 



been faster than in the early years of ACS, already caus- 
ing significant (~ 20%) signal loss in moderately bright 
sources of ~ 1000 e~/read, and ~ 50% for somewha t 
fainter sources of ~ 300 e"/read (jNoeske et al.l 12012ft . 
This faster degradation is believed to be due to the ex- 
treme solar activity minimum, and resulting cosmic ray 
maximum, during the initial flight years of UVIS. The 
resulting loss of data quality can be partially mitigated 
by post-processing. The effect is worse for very faint 
sources, whic h can be completely lost to "traps" be- 
fore r eadout (|MacKentv fc Smithl l2012t lAnderson et al.l 
120121) and cannot be recovered later. In the literature, 
CTE degradation is often referred to and mea sured as 
charg e transfer inefficiency (CTI = 1-CTE) (e.g. IMassevI 
2010), and we use this terminology interchangably below. 

After Epoch 2 of the UVUDF had already been 
obtained, the Space Telescope Science Institute 
STScI) released a new report on mitigating CTI 
MacKentv k, S mith 2012j). The strong recommendation 
is to use the "post-flash" capability of the instrument 
to illuminate the detector and add background light to 
the observation. This additional background will fill the 
traps and ensure that faint objects are not lost, as well 
as significantly improve the accuracy of pixel-based CTE 
corrections. This benefit comes at the cost of decreased 
sensitivity, however, due to the noise introduced by the 
added background. 

Considering that many of the science goals of the 
UVUDF rely on measuring (or setting limits on) the 
faintest sources, and require accurate photometry, we 
chose to follow the recommendation for post-flash. In 
Epoch 3, we applied a post-flash to bring the average 
background (the sum of post-flash, sky and dark cur- 
rent) up to about 13 electrons per pixel. In practice, 
this meant adding lie" in F225W and F275W, and 8e" 
in F336W. T he spatial distribution of post-flash light is 
not u niform ([MacKentv Sz S mith 20T2r lAnderson et al.l 
I20121 ). so target levels were set to ensure both a rea- 
sonable average and a sufficient background in the less 
illuminated regions. 

3.2. Binning the CCD Readout 

Without post-flash, the UVIS detectors are read-noise 
limited in the F225W and F275W filters, even in long ex- 
posures such as those needed for the UVUDF. The noise 
from the readout and from the sky background is about 
equal in F336W. As a consequence, there is the potential 
for tremendous sensitivity gain by binning the CCD pix- 
els 2 x 2 during readout. In principle, 2x2 binning results 
in a gain of a factor of 2 in signal to noise (S/N) ratio, 
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TABLE 2 
UVUDF Sensitivities 



Filter 


Zero Point a 


Epoch 


Exposure Time 


5a 0'.'2 ETC b 


5a 0'.'2 RMS 


50% completeness 




(mag) 




« 


(mag) 


(mag) 


(mag) 


F225W 


24.0403 


1&2 


39278 


28.3 


28.3 


28.6 


F275W 


24.1305 


1&2 


39106 


28.5 


28.4 


28.6 


F336W 


24.6682 


1&2 


45150 


29.0 


28.7 


28.9 


F225W 


24.0403 


3 


44072 


27.8 


27.9 


27.7 


F275W 


24.1305 


3 


41978 


27.7 


27.9 


27.7 


F336W 


24.6682 


3 


37646 


28.3 


28.3 


28.2 



Note. — UVUDF niters, zeropoints, and sensitivities. 
a Zeropoint information is available at http://www.stsci.edu/hst/wfc3/phot_zp_lbn 
b Exposure Time Calculator (ETC) modified to work with binned and post-flashed data. 



or 0.75 magnitudes. One concern in the decision to bin 
the CCD readout is the loss of spatial resolution. How- 
ever, the large number of repeated observations allow for 
excellent sub-pixel image reconstruction. 

Once the post-flash capability became available to mit- 
igate the effects of CTI, the benefit of binning the CCD 
readout was greatly reduced. At that point, the signal- 
to-noise gain would be under 20%, while still reducing 
the spatial resolution. As a result of these considera- 
tions, we chose to take the second half of the UVUDF 
data, that is Epoch 3, without binning the CCD readout. 

3.3. Parallel ACS Observations 

Coordinated parallel exposures were obtained with the 
ACS/WFC3 during the primary WFC3/UVIS observa- 
tions. Figure [3] shows the location of the parallel fields 
in comparison to other data in GOODS-South. Table [T] 
gives the specification for each parallel field, with posi- 
tion angle specified by the HST ORIENT keyword, which 
is the position angle of the U3 axis, where U3= — V3. 
The V3 angle is an angle based on the reference frame 
of the telescope, where V3 is perpendicular to the solar- 
array rotation axis. This angle describes the angle of 
rotation of the WFC3 UVIS exposure on the sky, and 
the position and rotation of the parallels. 

The Epoch 1 parallel exposures fall within the ERS 
field. The Epoch 2 parallel exposures fall outside 
of the main CANDELS and GOODS footprint, but 
still within the field observed by the GEMS program 
( Rix et al1l2004D. an d the ground-based U- and R-bands 
( Nonino et al.l l2009). Scheduling constraints did not per- 
mit a more favorable orientation. The Epoch 3 orienta- 
tion was chosen specifically to place the parallel field at 
the position of the parallels of the HUDF09. The distri- 
bution of exposures per filter varies with the position of 
the parallel data (see Table [TJ). 

The parallels of Epoch 1, which fall within the ERS, 
consist of 18 orbits. Given the depth of existing data 
in that field, we chose to obtain images with the three 
standard optical filters F435W, F606W, F814W. Four or- 
bit depth was obtained in the B-band (F435W), to more 
than double the previous imaging exposure time. The 
V (F606W) and wide Fband (F814W) exposure times 
were chosen so that when combined with previous imag- 
ing, the ratio would be ~1:2, following the strategy for 
parallel imaging in CANDELS. 

For Epochs 2 and 3, we obtained very deep B-band 
imaging. There is a growing recognition that HST's 
UV and blue optical capabilities are a unique resource 



which should be used now to prepare for later years 
when space-based observing will be focused on the near- 
infrared, w ith missions su c h as the James Webb Space 
Telescope (jGardner et all [2006h. the Wide Field In- 
frare d Survey Telescope flDressler et alj 12012). and Eu- 
clid (jLaureiis et all 120121 ) . With 20 and 46 orbits in 
Epochs 2 and 3 respectively, we obtained deep-field qual- 
ity images. For comparison, the original ACS/WFC3 B- 
band images of the UDF were obtained in 56 orbits. We 
note, however, that the detector performance was bet- 
ter at that time (see Sections 13.11 and I5.ip . In Epoch 2, 
we also obtained shallow imaging in the V (F606W), i 
(F775W), and z (F850LP) filters, to augment the shal- 
lower imaging available from GEMS. The failed visits de- 
scribed in section [3] shifted 4 orbits from planned B-band 
exposures in Epoch 2 to Epoch 3. 



HUDF09-1 



Parallel 1 



Parallel 2 
HUDF09 P1 
HUDF09-2 



HUDF05P34 




ERS IR 



Parallel 3 



CANDELS 
Deep IR 



CANDELS 
Wide IR 



GOODS South 

Fig. 3. — The footprint of the the ACS parallel pointings 
for Epochs 1, 2, and 3 are shown as purple squares. The 
greyscale image i s the V-band ACS map of GOODS-South from 
IGiavalisco et al.l 1120041 ), with North up and East left. The blue 
squares and nearby shaded regions indicate the footprint of the 
UVUDF UVIS pointings and complementary data from Figure [2] 
The blue shaded region indicates the footprint of the ERS imaging 
(Windhorst et al. 2011), the purple and brown indicate the foot- 
print for CANDELS Deep and Wide respectively HGrogin et al.l 
120111 1 , and the shaded red regions indicat e the footprint of th e 
near-infrared imaging from the HUDF09 (Bouwcns ct al. 2011). 
The green shaded region indicates the footprint of the HUDF09 
parallel 1, a nd the cyan shade d regions represents the HUDF05 
parallel P34 (Ocsch ct al. 2007). 
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4. DATA REDUCTION 

The UVUDF data set consists of four exposures and 2 
orbits per visit, with visits divided into three observing 
epochs as described above. In total, there are 15 visits 
(30 orbits) for each of the three filters. 

In this section, we describe the data reduction pro- 
cess needed to produce science quality images from the 
UVUDF observations. We plan to release fully reduced 
images and catalogs at a later time, but not in combi- 
nation with this paper (see Section [7]). Nonetheless, it 
is important to document the many issues with the data 
from this HST Treasury program, and for the reader to 
understand the process that led to the images used for 
the analysis in the later sections of the paper. The same 
lessons learned here will be relevant to planning of future 
UVIS observations. 

Binned and unbinned data (and data with and with- 
out post-flash) must be processed differently, and they 
require different calibration files. The software pipeline 
that we use for UVUDF data begins with the standard 
Pyraf/STSDA93 CALWF3 modules, though calibration 
files needed to be constructed with some care as de- 
scribed in this section. The processing of ACS parallel 
data closely follows the procedures used by CANDELS 
(|Koekemoer et aLll2011[) . and is not described here. 

4.1. Calibration Pipeline 

Calibration exposures (darks, biases, flat fields) for 
UVIS are obtained by STScI as part of the standard cal- 
ibration observations. In most cases, these calibrations 
are taken without binning the CCD readout, though a 
few binned calibration observations have been obtained. 
The CCD detectors are periodically heated in order to 
mitigate hot pixels that develop over time, called anneal- 
ing. Specifically, ~ 500 new hot pixels appear per day, 
while the annealing proces s removes > 70% of hot pixels 
([Borders fc Baggctt 20Q9|). The number of permanent 
hot pixels that can not be fixed by anneals is growing by 
0.05-1% per month (WFC3 instrument handbook). In 
order to minimize the number of hot pixels at any given 
time, the detector is annealed once per month. 

New calibration files are needed for the calibration 
pipeline: new biases, darks, and flats for the binned data, 
and new darks for the unbinned data. Only the bias 
files used data that were taken with onboard binning. In 
the other cases, unbinned calibration data are the basis 
of creating new files, with after-the-fact binning applied 
where necessary. We validated this latter procedure us- 
ing the limited set of onboard-binned calibrations that 
are available. We use a combination of custom scripts 
and standard STSDAS routines to make these calibra- 
tion files. The steps involved to construct each type of 
calibration file are described below. 

The standard calibration pipeline begins with an over- 
all bias correction, calculated by fitting the overscan 
region in a master bias frame, and removing the elec- 
tronic zero point bias level. Next, a bias reference frame 
is subtracted from the full image to correct for pixel- 
to-pixel bias structure. For the binned Epoch 1 and 
2 data, this reference file is created using the STScI 

18 Further documentation fo r all the PyRAF/STSDAS da ta re- 
duction software is provided at http://stsdas.stsci.edu/ 



software wfc3_reference.py ([Martel et al.l [20081) to aver- 
age ten onboard-binned bias frame exposures (Baggett, 
CAL-12798). For the unbinned Epoch 3 data, we use the 
standard unbinned bias frames provided by STScI. 

The next calibration step is the subtraction of a dark 
reference file to correct dark current structure and to 
mitigate hot pixels that can cause significant artifacts in 
the images. STScI releases new darks every 4 days that 
are based on the average of ~ 10 — 20 dark exposures 
with integration times of ~ 900s each. This is necessary 
due to the large number of new hot pixels per day, and 
the drastic change in hot pixels after each anneal. How- 
ever, binned darks are not obtained on a regular basis. 
Therefore, unbinned darks are binned after the fact us- 
ing custom IDL scripts. We validate this approach by 
measuring the dark current in one set of binned dark ex- 
posures obtained for this purpose (Baggett, CAL-12798). 

We find the standard processing of the dark calibra- 
tion is insufficient for the UVUDF data. The STScI- 
processed darks were created with two choices that are 
not optimized for this case. First, the process uses an 
unagressive definition of a hot pixel as a ~ 10c devia- 
tion. The choice results in warm-to-hot pixels not being 
masked in the UVUDF images, which add significant ar- 
tifacts to the highly sensitive mosaics. Second, the stan- 
dard processing uses the median value of the average 
darks (with hot pixels masked) as the value of all pixels 
in the dark frame. This median- value dark with hot pix- 
els is the calibration file available from STScI. It is not 
suitable for UVUDF data, because there is a low-level 
gradient present in the dark that is not subtracted. This 
gradient is typically small compared to the sky back- 
ground in the optical, with a peak-to-valley deviation of 
~ 3 e~/pix/hr. However, in the low-background NUV 
images, it is the dominant structure. While this back- 
ground can be corrected in the background subtraction 
phase of the pipeline, it is more accurate to model this 
background and subtract it before dividing by the flats. 

We therefore reprocess the darks, starting with the raw 
dark observations, using a procedure based on the one 
provided to us by STScI (J. Biretta, private communica - 
tion) and the wfc3_reference.py code (jMartel et al.l2 008). 
We ma ke two significant modifica tions to the STScI pro- 
cedure (jBorders fe Baggett 120091 ) to fix the issues iden- 
tified above. First, we use an iterative ~ 3er cutoff for 
defining a hot pixel, applied to cosmic-ray cleaned darks 
made from the average of a minimum of 10 exposures. 
This change significantly increases the number of hot pix- 
els masked (~ 7% of the image), but decreases the extra 
systematic noise. Secondly, we fit a 7 th order polyno- 
mial to the remaining non-flagged pixels in the image of 
each UVIS CCD. Then, we create a final dark frame by 
combining this polynomial fit, as the background value, 
and the hot pixels superimposed and flagged in the data 
quality array. In the case of the binned data, the new 
darks are binned after the fact. 

The next calibration step is the application of the flat 
field reference files. For the binned data, flats arc binned 
after the fact from the unbinned calibration data. For 
the unbinned data, the flats provided by STScI are ap- 
plied. The final calibration step is populating the pho- 
tometry keywords in the FITS header using the current 
filter throughput curves and detector sensitivity informa- 



tion using calwf3. 

The last processing step is the background subtraction 
of the individual calibrated images. The unbinned data 
have an artificial background introduced by the post- 
flash process. We subtract the post-flash reference files 
provided by STScI from the unbinned data. These ref- 
erence files are generated by STScI from stacks of post- 
flashed exposures, and then scaled to the flash count rate 
when applied to the data. However, both these post- 
flash-subtracted images, and the binned images, have 
a residual nonuniform background. We therefore fit a 
background to each individual image via a custom in- 
verse distance code. This code masks large fractions of 
each image for cosmic rays, sources, hot pixels, and bad 
pixels. It then interpolates the background value at any 
given pixel based on an inverse-distance weighting within 
a subgrid region. These backgrounds are then subtracted 
from all science images. The final products of the calibra- 
tion pipeline are basic calibrated background-subtracted 
images, together with data quality maps, that can be 
used as input to the mosaicking pipeline. 

Image registration and mosaicking are performed fol- 
lowing the procedures used for CANDELS. We refer 
the interested reader to Koekemoer et al. (2011). 
UV mosaics are regis tered to the ACS B-band image 
(jBeckwith et al J 12009) . 

4.2. Object Detection and Photometry 

We u se the Source Extractor s oftware version 2.5 (SEx- 
tractor; iBertin fe ArnoutHll996[ ) for object detection and 
photometry. SExtractor is used in dual image mode, 
where objects are detected in the deeper F435W (B- 
band) mosaic (Bcckwith et al. 2006), and the photom- 
etry is measured in a combined Epoch 1 and 2 mosaic 
and Epoch 3 mosaic for each filter. In this way, colors of 
sources are measured using the same isophotal apertures, 
and fluxes are measured for all B-band detected objects 
regardless of any flux decrement in the NUV mosaics due 
to the Lyman Break. Edge regions and the central chip 
gaps of the mosaics are excluded, and are set to the sky 
level with the same noise properties as the mosaics such 
that SExtractor does not find spurious sources along the 
edges or in the central chip gap. 

The detection parameters for the B-band mosaic are 
tuned such that no sources are detected in the nega- 
tive image. This is accomplished by setting the mini- 
mum area of adjoining pixels to 9 pixels, and a la de- 
tection threshold. A Gaussian filter is applied on the 
mosaics, with a full width at half maximum (FWHM) of 
3 pixels for object detection. SExtractor is provided an 
RMS weight map for both the detection and analysis im- 
age. The gain parameters are set to the exposure time, 
such that SExtractor calculates the uncertainties prop- 
erly. All source photometry has the the local background 
subtracted by SExtractor, using a local annulus that is 
24 pixels wide (with the inner radius depeding on source 
size). Zero points of 24.0403, 24.1305, and 24.6682 are 
applied for the F225W, F275W, and F336W mosaics re- 
spectively (see Table [2]). We note that since the B-band 
is significantly more sensitive than the NUV images, the 
resulting catalog contains B-band objects too faint to be 
measured in the UV, and thus cuts on the catalog are 
used as needed for each scientific purpose. 

The photometry of objects is measured with SEx- 



tractor using both isophotal and iKronl (|1980[ ) ellipti- 
cal apertures. Isophotal apertures are used whenever 
measuring the color of a source, such as in the color- 
color selections used in section 16.11 For this purpose, 
we also run SExtracto r on the F606W (V-band) mosaic 
(jBeckwith et al.ll2006D in dual image mode, still using the 
B-band as the detection image. This procedure results 
in aperture-matched photometry, although it is not cor- 
rected for variation in the point-spread function (PSF). 
Since the PSFs of the NUV and the optical B- and V- 
bands are quite similar, this correction will be small for 
these bands. For this overview paper, these color mea- 
surements are sufficient. Uncertainties will be dominated 
by CTI effects (see section [5~Tjl . Kron elliptical apertures 
are used to measure the total magnitude of each source 
via SExtractor's MAG_AUTO parameter. These magni- 
tudes measure the total flux from a source, and are used 
whenever a total magnitude is needed, such as in the 
number counts of LBGs (see section [53]) . 

5. DATA CHARACTERIZATION 

5.1. CTI effects 

Radiation damage sustained by the CCD degrades its 
ability to transfer electrons from one pixel to the next, 
trapping electrons (in part temporarily) during readout, 
while other electrons are moved to the next pixel. This 
results in trails of electrons in the direction of the CCD 
readout, with regions of the CCD furthest from the read- 
out affected most severely. The three different orienta- 
tions of the three UVUDF epochs enables the measure- 
ment of CTI effects in the data. Specifically, Epoch 1 and 
2 are at an angle of 101.25 degrees relative to each other, 
resulting in some galaxies located close to the readout 
in one epoch, and far from the readout in the other (see 
Figure^]). This configuration allows the characterization 
of the effect of CTI on the photometry and morphology, 
as well as an estimate of the number of faint galaxies that 
are completely lost. 

5.1.1. Corrections for CTI 

There are currently two methodologies to correct 
the photometry for CTI losses. The best method is 
a pixel-based CTE correction of the raw data based 
on empirical modeling of hot pixels in dark exposures 
(jAnderson fc Bedinl [20101 iMassev etaLl 120101) . Such a 
correction not only corrects the photometry, but also 
restores the morphology of sources (see section 15 . 1 ,3|) . 
A preliminary version of software to implement such a 
correction for unbinned WFC3/UVIS data was released 
in March 20lJ ra l. but significant improvements and ver- 
ification will be needed before the correction is stable 
enough to warrant the public release of corrected high 
level science products for the UVUDF. There are three 
major issues that need to be overcome for the software 
to fully support the UVUDF data: 1) The code only 
works for unbinned data, and half the UVUDF data are 
binned. 2) The current algorithm over predict s the CTE 
corre ction for low background faint sources ([Anderson! 
120131 ). and the binned half of the data have very low 
backgrounds. 3) Read noise mitigation in the algorithm 

19 For more information about the pixel- 
based CTE correction for WFC3/U VIS, see 
http : //www. stsci . edu/hst/wf c3/tools/cte_tools 
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resul ts in under-correction for faint sources (Anderson 
12013ft . The WFC3 team at STScI is aware of the lat- 
ter two limitations and is working on improvements. In 
addition, while the post-flash Epoch 3 data can have 
the CTE algorithm applied in a straightforward manner, 
post-flashed CTE corrected darks are required to match 
the hot pixels. 

The second method to mitigate the effect of CTI 
is to apply a correction to the measured flux densi- 
ties of sources, based on their location on the detector, 
the observation da t e, and their flux in electrons (e.g. 
Cawlev et all [2001] IRiess fc Mack! [200l IRhodes et all 
2007tlNoeske et all 120121 Bendregal et al. 2013). How- 
ever, the current WFC3 UVIS implementat ion of this 
catalog-based calibration (|Noeske et all [20121) has many 
limitations. First, it can only be applied for a small num- 
ber of quantized background levels, including virtually 
no background, ~ 3e _ /pix, and 20-30e~/pix. Thus, it is 
only applicable to the UVUDF Epoch 1 and 2 data for 
F275W and F225W, and these corrections have slightly 
higher backgrounds than the UVUDF data. The F336W 
data and all the Epoch 3 data have backgrounds that are 
not similar to any of the standard calibrations. The poor 
sampling of background levels in the calibrations makes 
interpolating between them unreliable. Second, the cali- 
bration was measured for relatively bright point sources, 
and the correction is uncertain at the faint end, which 
encompasses the majority of UVUDF sources. Third, it 
does not take into account other nearby sources which 
can fill charge traps and thereby shield the sources, 
lastly, it does not take into account the morphology 
of sources (i.e. size and shape), and therefore does 
not account for effects such as self-shielding that accom- 
pany non-point sources, as electrons from the part of the 
source closer to the readout will shield the other part 
from charge traps. 

K eeping these s everal limitations in mind, we apply 
the iNoeske et ail (|2012ft correction to the F225W and 
F275W mosaics of Epochs 1 and 2 separately. This cor- 
rection enables us to refine our investigation of the effects 
of CTI (e.g. Sections 15.1.21 and I5.1.4[) . However, we can 
not apply the calibration to the combined Epoch 1 and 
2 mosaic nor the F336W mosaics, so the science inves- 
tigations in Section 16.11 do not include the correction. 
Those investigations use the combined Epoch 1 and 2 
mosaic, which partially mitigates the CTI effect because 
objects far from the readout in one epoch are averaged 
with their counterparts closer to the readout in the other 
epoch. Future work using this data will apply the pixel- 
based CTE correction (when it is stable) to obtain more 
reliable photometry. 

5.1.2. CTI effects on photometry 

In order to characterize the CTI in the UVUDF data, 
a new catalog was created, with a method differing from 
that described in section |4~21 For each single epoch mo- 
saic in each filter, SExtractor was run in dual image 
mode, with the combined Epoch 1 and 2 mosaic as the 
detection image. The detection threshold was set such 
that we do not detect sources in the negative image. Ob- 
jects near the edges or near the chip gaps for any of the 
three epochs were excluded, and objects were required to 
be covered by all three epochs of observation. The cata- 
log was trimmed to only include sources with S/N ratios 



greater than 5er in all three single epoch mosaics. Galax- 
ies in the NUV images are often clumpy, which results 
in single galaxies appearing as multiple clumps in the 
images. Regardless of the deblending parameters used 
with SExtractor, these galaxies are detected as separate 
sources. This is not an issue for the CTI measurements 
described below, and the main catalog is not strongly af- 
fected by this, since the B-band is used as the detection 
image in that case. 

The effects of CTI are worst i n exp osures with low 
background ( MacKentv fc Smith! 12012ft , thus the mea- 
sured UVUDF CTI effects are described for F275W, 
which has a lower background than F336W, yet sources 
are brighter than in F225W, enabling us to measure more 
sources. Specifically, the unbinned equivalent average 
backgrounds are ~ 5.8 e~/pix/hr, ~ 6.2 e _ /pix/hr, and 
~ 12.2 e"/pix/hr for F225W, F275W, and F336W ex- 
posures, corresponding to ~ 2.4 e~/pix, ~ 2.5 e~/pix, 
and ~ 5.1 c~/pix in the half orbit exposures used. The 
backgrounds in F225W and F275W are consistent with 
the expected value due to dark current. The CTI effects 
are present in all three bands, but expected to be at 
a lower level in the higher-background F336W mosaics. 
The basic effect of CTI on the photometry is that the 
objects lose a fraction of their flux proportional to their 
distance from the readout, as electrons encounter more 
charge traps the further they travel. 

The uncorrected photometry of Epochs 1 and 2 are 
compared in the top panel of Figure [H which plots the 
difference in isophotal magnitude of Epoch 1 and 2 as a 
function of the Epoch 1 isophotal magnitude. The scatter 
is much larger than the expected la dispersion (shown as 
the gray shaded region) likely due to the effects of CTI. 
For objects far from the readout, the actual la dispersion 
is larger by more than a factor of 2 than expected. The 
photometric scatter is characterized as a function of the 
difference in distance to the readout between the epochs, 
as measured on the drizzled images. When the differ- 
ence is a large negative number, the sources are close 
to the readout in Epoch 1 and far from the readout in 
Epoch 2 (open blue squares). When the difference is a 
large positive number, the sources are close to the read- 
out in Epoch 2 and far from the readout in Epoch 1 (red 
filled circles). If CTI is the cause of the large scatter, 
then the expected behavior is for the blue squares to be 
primarily below the zero line, and the red circles to be 
primarily above the zero line. This behavior is indeed 
what is observed, confirming that CTI is the most likely 
cause of the large observed scatter. 

The CTE-corrected photometry of Epochs 1 and 2 are 
compared in the bottom panel of Figure [4l which plots 
the same quantities as the top panel, with the addition 
of a catalog-based CTE correction (see Section 15.1.11) . 
The CTE correction reduces the scatter observed in the 
top panel, and it removes the systematic offset of the 
red circles furthest from the readout. However, the scat- 
ter remains larger than expected, possibly due in part to 
the limitations of the catalog-based CTE corrections de- 
scribed in Section 15.1.11 On the other hand, the scatter 
could result from imperfect image registration, or CTI 
effects on source morphology causing inappropriate aper- 
tures to be used in the photometric measurements (see 
Section 15.1.31) . The image registration is unlikely to be 
the cause, because the Epoch 1 and 2 have relative as- 
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trometric accuracy of better than 0'.'05. It is possible 
that the CTI effects on source morphology is the cause, 
though the use of the combined Epoch 1 and 2 mosaic 
as the detection image somewhat reduces this effect (but 
see Section 15 .1.3[) . 

We test the hypothesis that something other than CTI 
is the cause of the scatter by making a comparison that 
is mostly insensitive to the distance to the readout. We 
compared two subsets of the Epoch 2 data, each consist- 
ing of half the exposures (2a and 2b). Figure [S] plots 
the difference in isophotal magnitude between the two 
halves of the Epoch 2 data as a function of the Epoch 
2a isophotal magnitude. The points are color coded 
by distance to the readout, as no difference in read- 
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Fig. 4. — Photometry comparison of sources in Epochs 1 and 
2 F275W mosaics illustrating a larger photometric scatter than 
expected from the uncertainties, likely due to CTI effects. For ob- 
jects far from the readout, the actual lcr dispersion is larger by 
more than a factor of 2 than expected. Both panels are the same 
except that the bottom panel includes a catalog based CTE cor- 
rection for both E pochs 1 and 2 assuming point source morphology 
(see Section [5. 1,111 , The difference in isophotal magnitudes between 
Epochs 1 and 2 should be zero, with a scatter that increases with 
increasing magnitude. The black line is the zero difference line, and 
the expected lcr dispersion is shown as the gray shaded region from 
uncertainties as measured by Source Extractor. The colors denote 
the difference in source distance to readout between the epochs. 
The blue open squares are sources close to the readout in Epoch 1 
and far from the readout in Epoch 2, the green filled triangles are 
sources an intermediate distance from the readouts, and the red 
filled circles are sources far from the readout in Epoch 1 and close 
to the readout in Epoch 2. 
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Fig. 5. — Photometry comparison of Epochs 2a and 2b F275W 
with a photometric scatter mostly consistent with that expected 
from the uncertainties. No correction for CTI is applied, because 
the correction would be a function of distance to the readout and 
therefore the same in both halves of the epoch 2 data. The expected 
lcr dispersion is shown as the gray shaded region. The blue open 
squares are sources close to the readout (< 700 pixels away), the 
green filled triangles are sources an intermediate distance from the 
readout (> 700 pixels and < 2200 pixels away), and the red filled 
circles are sources far from the readout (> 2200 pixels away). Re- 
gardless of distance to the readout, source magnitudes are mostly 
consistent with the lcr scatter (gray shaded region). 



out distance exists. Regardless of the distance to the 
readout, magnitude differences are consistent with ran- 
dom scatter, although with a slightly larger magnitude 
than expected from the measurement uncertainties (gray 
shaded la dispersion). This minor remaining difference 
is most likely due to a slight underestimation of the un- 
certainties by SExtractor, possibly caused by SExtrac- 
tor not including the uncertainty in local sky subtrac- 
tion. It has been noted several times in the literature 
that SExtractor underestimates the true uncertainties 
(iFeldmeier et alj|2002t lLabbe et al.ll2003t iGa wiser et al.l 
l2006HBecker et alJl2007t iCoe et al.M2013f ). 

Another method to visualize the CTI effects is to plot 
the magnitude difference in Epochs 1 and 2 versus the 
difference in distance to the readout (top panel, Figure 
[5]). Sources falling to the left in this figure are close to the 
readout in Epoch 1 and far from the readout in Epoch 2, 
while sources falling to the right in this figure are close 
to the readout in Epoch 2 and far from the readout in 
Epoch 1. Sources for which electrons travel larger dis- 
tances to the readout lose more flux, so CTI effects would 
cause the difference in magnitude to be negative in the 
left half of the figure and positive in the right side of the 
figure. The sources used in the figure are color coded by 
magnitude, with purple triangles representing the bright- 
est, and green circles representing the faintest. The pur- 
ple triangles have a smaller scatter, consistent with the 
fact that bright sou rces are less s everely affected by CTI 
than faint sources (Masscy 2010). The red points with 
error bars in Figure [6j which show the average values in 
bins of equal numbers per bin, emphasize the trend. The 
uncertainties are the standard deviation of the points in 
each bin divided by the square root of the number of 
points per bin. The bottom panel of Figure [5] is the 
same as the top, with the addition of a catalog-based 
CTE correction (see Section [5. 1.1[) . The CTE correction 
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Fig. 6. — Photometry comparison of sources in Epochs 1 and 
2 F275W mosaics as a function of the difference in distance to 
readout. Sources falling to the left in this figure are close to the 
readout in Epoch 1, and far from the readout in Epoch 2, while 
sources falling to the right in this figure are close to the readout 
in Epoch 2, and far from the readout in Epoch 1. Sources are 
color coded by their Epoch 1 magnitudes. The black line is the 
zero difference line. The red points with error bars are average 
binned values, with equal numbers of sources in each bin. Observed 
photometry is consistent with CTI effects, with the difference in 
magnitudes being negative in the top side of the figure, and positive 
in the right side of the figure. The slope of the effect is removed 
(bottom pa nel) w hen applying the catalog based CTE correction 
(see Section |5.1.1| I. 

somewhat reduces the scatter observed in the top panel, 
and removes the slope observed in the data. 

Given that the increased photometric scatter is corre- 
lated with the readout direction, and that there is signif- 
icantly less scatter when comparing the subsets of Epoch 
2 data, we conclude that CTI is the dominant cause of 
the large scatter in photometry observed in Figures |4] 
and [3 It is possible that other calibration issues con- 
tribute as well, but they would require effects that are 
also dependent on source position on the detector. 

5.1.3. CTI effects on morphology 

CTI affects the shape of g a laxy i mages as well as their 
photometry. Rhodes e t al.l (|2010f ) investigated the ef- 
fects of CTI on galaxy morphology using simulations and 
found that small galaxies are more affected by CTI than 
large ones. They also found that small bright galaxies are 



slightly less affected by CTI than small faint ones, but 
this dependence is not observed for large galaxies. The 
net effect of CTI on image morphology is a smearing out 
of the flux in the readout direction. Thus CTI results 
in circular objects appearing elongated in the readout 
direction. 

This elongation effect is observed in the UVUDF data, 
as shown in the example in Figure [7] This galaxy is 
located about two thirds the length of the detector away 
from the readout in Epoch 1, and almost as far as possible 
from the readout in Epoch 2. In this example, both 
the bright galaxy and the nearby smaller structures are 
elongated in the readout direction, as marked by the red 
lines. The near 90 degree separation of Epochs 1 and 2 
shows the magnitude of the elongation in each direction. 
We note that this elongation also affects the astrometry, 
limiting the precision of the alignment between epochs 
and between UVIS and ACS images. 



Epoch 2 




B-band 



Fig. 7. — Example of a galaxy affected by CTI in the F275W 
images. The top left panel is from Epoch 1, the top right panel 
is from Epoch 2, the bottom left panel is from Epoch 3, and the 
bottom right panel is the B-band image. The red arrows corre- 
spond to the readout direction, and the galaxy is elongated in the 
readout direction in each case. The elongation is worst in Epoch 
2, as it is furthest from the readout in that Epoch. The elongation 
is reduced in the post-flashed Epoch 3 compared to Epoch 1. 



5.1.4. CTI effects resulting in source loss 

Another effect of CTI on the images is the possibil- 
ity of losing faint sources completely. Studies of warm 
pixels in long dark exposures show that the number of 
warm pixels decrease drastically further away from the 
readout, and the effect is worse for fainter warm pixels 
(|MacKentv k, S mith 20121) . That study is a worst case 
scenario, because warm pixels are not shielded by other 
nearby pixels as is the case for pixels associated with faint 
astronomical sources. Nonetheless, post-flash calibration 
observations of Omega Centuri confirm that faint sources 
in low backgrounds ca n disappear completely due to CTI 
(|Anderson et al.ll2012l ). The sensitivity limit of observa- 
tions is thus set by the exposure time of each individual 
exposure rather than the average of a stack. This depth 
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Fig. 8. — Histograms of sources detected in the F275W Epoch 
1, 2, and 3 mosaics based on their isophotal CTE corrected mag- 
nitudes. Sources are split into two groups based on distance to the 
readout, with sources in the halves of the chips close to the read- 
out shown in blue, and the sources in the other halves far from the 
readou t sho wn in red. The 50% and 10% completeness levels (see 
section I5.2|l are plotted in green and brown respectively. Sources 
close to the readout appear to have a tail beyond the 10% com- 
pleteness, while sources far from the readout drop more steeply. 
This suggests that we are losing sources far from the readout that 
are not lost close to the readout. 



varies as a function of distance to the readout, morphol- 
ogy, and position of other sources on the detector. 
A simple test of source losses is a comparison of the 



number counts of detected objects as a function of mag- 
nitude for sources close and far from the readout. We 
start with the B-band selected catalog described in Sec- 
tion 14.21 and consider sources down to lOrr detections 
in the B-band and 3cr detections in F275W. Except for 
Epoch 3, we apply the catalog based CTE correction (see 
Sections [5.1.11 and 15 . 1 .2\ to reduce the effects of CTI on 
the photometry. The prediction is that some sources far 
from the readout in Epochs 1 and 2 will be lost com- 
pletely, while source losses should be greatly reduced in 
the post-flashed Epoch 3 data. 

A histogram of detected sources based on their isopho- 
tal CTE-corrected magnitudes is shown in Figure [5] for 
all three epochs. For Epochs 1 and 2, more faint sources 
are found close to the readout (blue) than are found far 
from the readout (red) , suggesting that some sources far 
from the readout have been lost. There is no significant 
difference in the number of B-band sources in the same 
sample areas. 

The sources lost due to CTI are very faint, and the 
number counts at these faint magnitudes are suppressed 
by the incompleteness due to lack of sensitivity (Section 
15. 2[) . It is therefore difficult to distinguish sources lost 
due to CTI from sources that would not be detected be- 
cause of insufficient sensitivity. Most of the losses from 
CTI are at magnitudes close to the 10% completeness 
limit for the Epoch 1 and 2 F275W mosaics. Keeping this 
limitation in mind, as well as the small number statistics 
at the faintest magnitudes where incompleteness is very 
high, we estimate the number of lost sources by compar- 
ing the number counts close to and far from the readout. 
For sources fainter than the 50% completeness limit of 
AB^28.3 mag (see section I5.2|) . we find that at least 
~ 30 sources are lost in each Epoch 1 and Epoch 2 (out 
of ~600 source positions that are common between the 
epochs), while no sources are lost in Epoch 3 (out of 
~ 500). The total number of lost galaxies is likely larger 
than those found above, because sources in the middle of 
the CCDs may also be lost. These sources are not close 
to the readout in either Epochs 1 and 2 and fall below 
the sensitivity limit of Epoch 3, making them difficult to 
identify. We expect that the number of these sources per 
area is smaller than the number found far from the read- 
out, suggesting the total number of lost sources is likely 
within a factor of two of those observed to be lost. Our 
best estimate is a loss of < 100 sources out of ^600. The 
small number of losses suggests that the results presented 
in Section 16.11 are not strongly biased due to CTI. 

Another empirical test of source losses due to CTI is 
to compare indvidual sources that are detected close to 
the readout in one epoch but whose position is far from 
the readout in another epoch. The Epoch 2 mosaic is 
slightly more sensitive than the Epoch 1 mosaic (8 orbits 
F275W in Epoch 2 compared to 6 orbits for Epoch 1), so 
sources that are detected in Epoch 1 close to the readout, 
but not detected in Epoch 2 far from the readout demon- 
strate the effect. In searching for such sources, we also 
required them to be det ected in the significa ntly more 
sensitive B-band image (|Beckwith et alJ [2006) . There 
exist a few such sources, and an example is shown in 
Figured The left panel is a cutout of the F275W Epoch 
1 mosaic and the middle panel is from Epoch 2, and the 
right panel is from Epoch 3. This source is observed in 
the optical ACS images, and is object 4188 in the catalog 
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bv iCoe et all (|2006ft . It has an F275W isophotal mag- 
nitude of 2 8.6 ±0.1, and a F435W total magnitude of 
27.9 ± 0.06 (|Coe et al.l l2006). This source is detected at 
8a in Epoch 1, and should have been observed at least 
at that S/N in the Epoch 2 data. 

The potential loss of faint objects is one of the primary 
reasons that we decided to use the post-flash option in 
Epoch 3. The other motivations for the post- flash include 
reducing other CTI effects and significantly improving 
pixel-based C TE corrections by taking data with higher 
backgrounds (M acKentv fe Smithll2012l ). 

The evidence that no sources have been lost in Epoch 3 
is encouraging, though the sensitivity limit is necessarily 
worse. While the F275W exposure time in Epoch 3 is 
about double that of Epochs 1 and 2 individually, Epoch 
3 is significantly less sensitive (see Section IBTS"]) . In fact, 
most of the sources that appear to be lost in Epochs 1 
and 2 due to CTI would not have been detected in the 
Epoch 3 mosaic in the first place. Thus there are very few 
examples of sources that were lost in the data without 
post-flash but are present in Epoch 3. One such example 
is shown in Figure 1101 This galaxy is observed in the 
optical ACS ima ges, and is object 8020 in the catalog by 
ICoe et all (|2006f ). It has an F275W isophotal magnitude 
of 27.8±0. 1, and a F435W total magnitude of 27.9 ±0.04 
(|Coe et al.M2006h . The source is detected at 9a in Epoch 
3, and would have been easily detected in both Epochs 
1 and 2 were it not for CTI. 

The CTI effects create a dichotomy between the first 
two UVUDF epochs and the third epoch. The combined 
Epoch 1 and 2 mosaic is more sensitive than the Epoch 3 
data, but suffers more from CTI, and some objects may 
be lost completely. Once pixel-based CTE corrections 
are applied, the Epoch 3 data will be the best charac- 




FlG. 9. — Example of a galaxy lost due to CTI in the F275W 
Epoch 2 mosaic. The left panel is a cutout of the F275W Epoch 
1 mosaic, the middle panel is from the Epoch 2 mosaic, and the 
right panel is from the Epoch 3 mosaic. The galaxy is present in 
the shallower Epoch 1 data which is close to the readout, and is 
not detected in the slightly more sensitive Epoch 2 data which is 
far from the readout. This galaxy is observe d in the optical A CS 
images, and is object 4188 in the catalog by Coc et al. ( 20061), It 
has an F275W isophotal magnit ude of 28.6 ± 0. 1, and a F435W 
total magnitude of 27.9 ± 0.04 JCoT ct al. 2006]). The galaxy is 
detected at 8<r in Epoch 1, and should have been observed at a 
higher significance in Epoch 2 were it not for CTI. The galaxy is 
also detected in the shallower post-flashed Epoch 3. 



terized NUV mosaic available. We agree with the STScI 
recommendation that future WFC3 UVIS observations 
that require very sensitive measurements use the post- 
flash. 




Fig. 10. — Example of a galaxy lost due to CTI in the F275W 
Epoch 1 and 2 mosaics, but preserved in Epoch 3 due to post-flash. 
The left panel is a cutout of the F275W Epoch 1 mosaic, the middle 
panel is from the Epoch 2 mosaic, and the right panel is from the 
Epoch 3 mosaic. The galaxy is present in the shallower Epoch 3 
data, and is not detected in the more sensitive Epoch 1 and 2 data. 
The galaxy is approximately in the middle of the chip in all three 
epochs. This galaxy is observ ed in the optical ACS images, and is 
object 8020 in the catalog bv lCoe et all (|2006|). It has an F275W 
isophotal mag nitude of 27.8 ± .1, and a F435W total magnitude 
of 27.9 ± 0.04 l lCoe et al.112009) . The galaxy is detected at 9cr in 
Epoch 3, and would have been easily detected in both Epochs 1 
and 2 were it not for CTI. 



5.2. Sensitivity 

We use two common methods to characterize the sen- 
sitivity of the UVUDF data. First, we measure the sky 
fluctuations of the images. Secondly, we measure the 50% 
completeness limit, measured by recovery of simulated 
sources placed in the science mosaics. The complete- 
ness test takes into account both the sky surface bright- 
ness and the spatial resolution of the mosaic s, yielding a 
good sense of the usable depth of an image (|Chen et all 
2002t ISawicki fc Thompson] l2005t iRafelski et all 120091: 
Windhorst et alll2011l ) 

The sky noise of each image is measured via the pixel- 
to-pixel rms fluctuations. These fluctuations are mea- 
sured in 51x51 pixel boxes at 1000 semi-random loca- 
tions, such that the boxes are entirely on the image, do 
not fall on a detected object, and the boxes do not over- 
lap other boxes. This technique is deisnged to be less 
sensitive to any residual gradient in the image than sim- 
ply using the rms of all unmasked pixels. The rms in each 
mosaic is the iterative sigma clipped mean of the rms in 
each box, which is determined with an iterative sigma 
clipped standard deviation. This rms is then multiplied 
by the noise correlation ratio to account for the corre- 
lated noise from drizzling the mosaics. The approximate 
correlation ratio of the UVUDF data is ~ 2.5 and ~ 1.5 
for the binned and unbinned data, r espec tively, based 
on equation 9 from Fruchter & Hookj (|2002f h These rms 
values corrected for correlated noise match the expected 
values from the rms images. The resulting 5a rms magni- 
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tudes (assuming 0'.'2 radius aperture) for the mosaics are 
tabulated in Table [2] These values are within 0.1 — 0.2 
mags of the 0'.'2 aperture, 5er magnitudes predicted by 
the STScI exposure time calculator modified for binning 
or post- flash (see Table [5]). 

We performed a standard completeness test to confirm 
the noise characteristics of the data by planting and re- 
covering simulated objects. This test does not take into 
account the loss of sources at the faint end due to CTI, 
and so the results of the test are an upper limit on the 
completeness. Specifically, the 50% completeness magni- 
tude limit due to noise is measured by planting Gaussian 
PSFs for a range of magnitudes in the mosaics at semi- 
random locations, and counting the fraction of sources 
that are recovered with SExtractor. The PSF FWHMs 
are matched to those measured in the data for each filter. 
Unresolved sources are sel ected from the pub lished cata- 
logs of stars in the UDF (jPirzkal et al.l 120051 ) . However, 
there are only a small number of identified sources bright 
enough in the NUV to be used for PSF determination. 
Three sources are used for F336W, and two sources are 
used for F225W and F275W. The sources are each reg- 
istered to their subpixel centers, normalized by the peak 
value, and coadded with a mean. The resul t ing P SFs are 
worse than measured bv lWindhorst et al.l (|201 If ) in the 
ERS, because half the data are binned and CTI affects 
the source morphology. For the combined Epoch 1 and 
2 mosaic, we measure PSF FWHMs of 0'.'133,0'.'133, and 
07127 for F225W, F275W, and F336W. For the Epoch 3 
mosaic, we measure PSF FWHM's of 07134, and 07121 
for F275W and F336W, and use the F275W PSF for the 
F225W PSF as it is not well determined. The locations 
of the planted sources are constrained such that they do 
not fall off the edges, fall on a real detected object, or 
fall on any previously planted source. The detection effi- 
ciency as a function of magnitude is shown in Figure 111! 
and the 50% completeness magnitudes are tabulated in 
Tabled 

6. INITIAL RESULTS 

In this section, we briefly present initial results, repre- 
sentative of those that will be possible with the UVUDF. 
We describe the color selection of galaxies using the Ly- 
man break technique, and we demonstrate the utility of 
the deep NUV images for morphological analysis. In both 
cases, these results are presented based on the combined 
Epoch 1 and 2 mosaics, without the application of the 
pixel-based CTE correction. We anticipate that future 
papers will improve the analysis once that correction is 
stable and can be confidently applied. 

6.1. Lyman Break Galaxies 

The selection of high redshift galaxies by the iden- 
tification of the strong Lyman break feature in their 
SED using broad-band photometry has b ee n ex- 
trem e ly successful (e.g. iSteidel et all ll996bTlal 1 1999 . 
" 2003 lAdelberger et al.l 120041: [Bouwens et al.ll200l 2006 . 
20Tq I201H: [Rafelski etafl l200a IReddv et alT 12008 ; 
Reddv fc Steidell 120091: IReddv et all I2012bl ). Although 



less precise than a full SED-fit such as those used in pho- 
tometric redshift estimates, the Lyman break identifica- 
tion is a standard in the literature. Here, we have taken 
a first look at selecting LBGs in the UVUDF at redshifts 
where the Lyman break falls in the NUV filters: 1.7, 2.1, 
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Fig. 11. — Detection efficiency of the Epochs 1 & 2 (top) and 
Epoch 3 (bottom) mosaics as a function of total magnitude. These 
are the recovery percentages of simulated point sources in the im- 
ages. The limiting magnitude is defined as the magnitude at which 
50% of the sources are recovered. The limiting magnitudes for 
Epoch 1&2 are 28.6, 28.6, and 28.9 mag, and for Epoch 3 are 27.7, 
27.7, and 28.2 mag, for F225W, F275W, and F336W (see Tabled. 
This completeness test does not account for source losses due to 
CTI. 



and 2.7 in F225W, F275W, and F336W, respectively. We 
directly compare th ese initial results with published re- 
sults from the ERS (Windhorst et al.ll2011| ). which used 
the same filters in shallower data (AB=26.9) over a larger 
area (about 50 square acrminutes). 

We implemen t the d ropout criteria used on ERS data 
bv lHathi et al.l (|2010h (hereafter H10) and lOesch et al.l 
(|2010aO (hereafter O10). These consist of both color- 
color criteria as well as S /N criteria for candidates to be 
considered dropouts. 

Faint stars were removed from the sample by posi- 
tion matching sources in the UVUDF catalog with the 
published catalog o f unresolved sources in the UDF 
(Pirzkal et~a^[2005T ). 25 sources are found to match 
this catalog (0.1" matching radius); 22 of these sources 
are identified as stars according to the criteria of 
iPirzkal et all (|2005f ). There are 7, 10, 18 stars detected 
at S/N = 3 threshold in F225W, F275W, F336W, re- 
spectively. 
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TABLE 3 

Summary of Dropout Galaxies in ERS and UVUDF. 



Dropout ERS UVUDF (ERS Depth) UVUDF (Full Depth) 







Predicted 


Observed 


Predicted 


Observed 


Predicted 


Observed 


Filter 


Method 


Number 


Number 


Density 


Number 


Number 


Density 


Number 


Number 


Density 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


F336W 


H10 


394 


256 


5.1 


49 


37 ±6.6 


6.0 ± 1.1 


185 


211 ± 14.5 


34 ± 2.3 


z ~ 2.7 


O10 


448 


403 


8.6 


56 


67 ±8.2 


10 ±1.2 


224 


304 ± 17.4 


49 ± 2.8 


F275W 


H10 


228 


151 


3.0 


28 


22 ± 5.2 


3.5 ±0.8 


86 


88 ± 9.4 


14 ± 1.5 


z ~ 2.1 


O10 


102 


99 


2.1 


13 


10 ±3.7 


1.6 ±0.6 


125 


146 ± 12.1 


24 ±2.0 


F225W 


H10 


62 


66 


1.3 


8 


4 ±2.5 


0.6 ±0.4 


36 


25 ± 5.5 


4.0 ±0.9 


z ~ 1.7 


O10 


99 


60 


1.3 


12 


9 ± 3.5 


1.5 ±0.6 


111 


61 ± 7.8 


9.8 ± 1.3 



Note. — Column (1) indicates the dropout filter and redshift bin. Col umn (2) indicates the re ference to the dropo ut method 
and luminosity function used to identify and predict source counts: H10 (Hathi et al. 2010]); O10 JOesch ct al. 2010a). Columns 
3-11) compare predict ed and observed source counts for each dropout type. ERS refers to the Early Release Science data 
Windhorst eFaUl2ljTlI) . UVUDF(ERS) refers to UVUDF data analyzed to a comparable depth as ERS, UVUDF(Full) refers 
to UVUDF data analyzed to its full depth. Dropout sky density values are in units of arcmin -2 . Uncertainties on the observed 
number and density of sources are Poissonian. Errors to predicted source counts are discussed in the text. 

6.1.1. F336W, F275W, F225W Dropouts 



For reference, the H10 criteria for dropout galaxies are 
given below. F336W dropouts require observed magni- 
tudes and signal-to- noise ratios, S/N, satisfy each of the 
following 



' F336W - F435W > 0.8 
F4351U < 26.5 
F4351U - F606W < 1.2 
F4351U - F606W > -0.2 
< F336W - FA35W > 0.35 + [1.3 x (FA35W - F606W)} 
S/N(F435W) > 3 
S/N(F336W) < 3 
S/N{F275W) < 1 
S/N{F225W) < 1 

(1) 

Similarly, F275W dropouts are identified by the crite- 
ria: 



' F275W - F336W > 1.0 
F336W < 26.5 
F336W - FA35W < 1.2 
F336W - FA35W > -0.2 

F275W - F336W > 0.35 + [1.3 x (F336W - F435W)] 
S/N(F336W) > 3 
S/N{F275W) < 3 
S/N{F225W) < 1 

(2) 

F225W dropouts require all of the following criteria: 



'F225W-F275W > 1.3 
F275W < 26.5 
F275W - F336W < 1.2 
F275W - F336W > -0.2 
' F225W - F275W > 0.35 + [1.3 x (F275W - F336W)} 
F336W - F435W > -0.5 
S/N(F275W) > 3 
S/N(F225W) < 3 

(3) 

Figures H21 H31 andlMlillustrate the dropout candidates 
selected according to the H10 criteria in color-color dia- 
grams. Sources with S/N < 1 in the dropout band have 
their fluxes replaced with la upper limits and are indi- 
cated as arrows in the figures. Stars are indicated as blue 
asterisks. Dropout candidates (defined as meeting all of 
the criteria of H10) are indicated as red symbols. 

For a direct comparison with previous observations, 
we use the UVUDF data to examine source counts at 
the sensitivity of the shallower ERS data and reproduce 
the selection at the shallow sensitivity. Because the H10 
and O10 dropout selection include a cut on source signif- 
icance, we scale the S/N ratio of the UVUDF measure- 
ments to what we would expect from the ERS, using the 
HST Exposure Time Calculator (ETC). The S/N ratios 
change by factors of 0.19, 0.33, 0.33 for F336W, F275W, 
F225W respectively. Dropout selection at the shallower 
limit is also shown in the figures. 

Table [3] indicates the number of dropout sources found 
using the methods of H10 and O10 in UVUDF data and, 
for comparison, those reported previously in ERS data. 
We find 37, 22, 4 H10 dropouts in F336W, F275W, 
F225W bands, respectively, to ERS depth in UVUDF, 
and we find 211, 88, 25 H10 dropouts in F336W, F275W, 
F225W bands, respectively, to the full depth in UVUDF. 
The raw number of dropouts in the narrow, deep UVUDF 
data is comparable to the numbers found in the wider, 
shallower ERS data. Table [3] also compares the sky den- 
sities of dropout candidates reported in H10 and O10 
to those detected in UVUDF. We find dropout sources 
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Fig. 12. — (left panel) Color-color diagram illustrating F336W dropout candidates in UVUDF to ERS depth using the method of Hathi 
et al 2010 (H10). Colors are computed from subtracting magnitudes in the appropriate HST wavebands. Sources are illustrated to ERS 
depth (F435W = 26.5 AB) and have S/N degraded to match ERS observations (as discussed in the text). Red circles are dropout galaxies 
according to criteria of H10. Stars are indicated as blue asterisks, and upper arrows indicate non-detections in the dropout band replaced 
with a la upper limit. Gray points indicate all sources in the UVUDF catalog to this depth limit; gray points in the color selection region 
fail to meet the S/N criteria of bona fide dropout galaxies. Solid lines indicate the color-selection region of H10; dashed lines indicate the 
color-selection region of O10. (right panel) Color-color diagram illustrating F336W dropout candidates to the full depth of UVUDF. Colors 
and symbols are the same as in the left panel. 




Fig. 13. — (left panel) Color-color diagram illustrating F275W dropout candidates in UVUDF to ERS depth using the method of Hathi 
et al 2010 (H10). Colors are computed from subtracting magnitudes in the appropriate HST wavebands. Sources are illustrated to ERS 
depth (F336 W = 26.5 AB) and have S/N degraded to match ERS observations (as discussed in the text). Colors and symbols are the same 
as in Figure [12] (right panel) Color-color diagram illustrating F275W dropout candidates to the full depth of UVUDF. 
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Fig. 14. — (Top panel) Color-color diagram illustrating F225W 
dropout candidates in UVUDF to ERS depth using the method 
of Hathi et al 2010 (H10). Colors are computed from subtracting 
magnitudes in the appropriate HST wavebands. Sources are illus- 
trated to ERS depth (F275W = 26.5 AB) and have S/N degraded 
to match ERS observations (as di scus sed in the text). Colors and 
symbols are the same as in Figure Il2l (Bottom panel) Color-color 
diagram illustrating F225W dropout candidates to the full depth 
of UVUDF. Colors and symbols are the same as in the top panel. 
There are 25 F225W dropouts. 



with comparable sky densities as H10 and O10 at the 
same depth and S/N limits. 

Table [3] also compares the dropout selection meth- 
ods of H10 and O10 applied in the UVUDF to each 
other; overall, the method of H10 is more conservative 
than the method of O10. For example, the number of 
O10 F336W dropouts and their resulting sky density ex- 
ceeds the number of H10 F336W dropouts by a factor 



of 1.8 (67 O10 dropouts vs. 37 H10 dropouts). This 
disparity arises for several reasons. First, H10 uses a 
F435W-selected catalog, whereas O10 uses a F606W- 
selected catalog; therefore, O10 dropout selection begins 
with a larger sample of sources (531 F606W detected 
sources for O10 versus 391 F435W detected sources 
for H10). Second, H10 has a more stringent require- 
ment for S/N in the bands blue- ward of the dropout 
band than O10 (S/N(F275W,225W) < 1 for H10 vs. 
S/N(F275W,F225W) < 2 for O10). Additional differ- 
ences between the two methods include the upper S/N 
limit for sources in the dropout band adopted by H10 
(S/N(F336W) < 3 for H10 vs. no S/N(F336W) require- 
ment for O10), the higher S/N requirement in the band 
redward of the dropout band in O10 (S/N(F435W) > 5 
for O10 vs. S/N(F435W) > 3 for H10), and the differ- 
ences in color selection regions between the two methods, 
as illustrated in Figure [T2J Stars are found and rejected 
in each sample with approximately the same percentage 
(15% and 20% for H10 and O10 respectively). 

6.1.2. Source Count Prediction 

We use the published luminosity functions of H10 and 
O10 dropout sources to predict an expectation for the 
number of sources to be found in the UVUDF. Each lu- 
minosity function, expressed as a space density of galax- 
ies, <f>, in units of Mpc~ 3 mag" 1 as a function of absolute 
magnitude, M, at rest-frame 1500 A is described by a 
Schechter function (|Schechterj|1976D . We use the fitted 
parameters M* and a reported in H10 and O10 for 
each dropout filter (F336W, F275W, F225W) to predict 
the space densities of galaxies in each redshift range. 

The differential number of galaxies per unit redshift, 
dN/dz, for each dropout filter, is given by integrals over 
the Schechter luminosity function, 3>(M), expressed as 
a function of absolute magnitude, M , multiplied by the 
(published in H10 and O10) gaussian galaxy redshift dis- 
tribution, g(z), and by the comoving volume element, 
dV/dzdfl, and finally integrated over the survey solid 
angle, ft. 



dN 

dz 



dQ 



dV 
dzdfl 



§{M)g{z)dM (4) 



The lower limit of integration is chosen to include the 
brightest observable galaxies. The number of sources 
down to the magnitude limit, Mn m , is found by inte- 
gration over the mean, z m , of the galaxy redshift distri- 
bution within ±lcr limits 



,+a dN 

N{< M Um ) = I —dz 



(5) 



No correction is made for completeness or selection ef- 
ficiency effects (i.e. an effective volume correction), since 
these corrections are specific to each H10 and O10 data 
set, and are not transferrable to the UVUDF data. 

The H10 and O10 luminosity functions were computed 
at rest-frame 1500 A, which for F225W, F275W, F336W 
dropout galaxies at redshifts z ~ 1.7, 2.1, 2.7, correspond 
to observed-frame 4050, 4650, 5550 A respectively. How- 
ever, H10, O10 and the UVUDF catalogs are selected 
from different wavebands (e.g. UVUDF uses a F4351V- 
selected catalog). To calculate the number density of 
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LBGs, the upper magnitude limit for ERS and UVUDF 
were modified by color-correction terms. These terms ac- 
count for the difference between rest-frame 1500 A and 
the catalog detection band. An estimate of the upper 
magnitude limit at rest-frame 1500 A is found by inter- 
polating between the limits for the two closest photo- 
metric bands. These correction factors, dm, were added 
to the catalog detection limits, mumn, to determine up- 
per limits of integration for the luminosity functions, 
miimit = mdetect+dm. For UVUDF data, we used a limit 
of rridetect — 28 to avoid the magnitude range in which 
sources can be lost to CTI. Correction factors were found 
to be dm = +0.130, -0.066, -0.236 for F225W, F275W, 
F336W dropouts respectively. 

For verification, we use the reported LF to predict the 
number counts that were observed in the ERS data them- 
selves, the expected number counts in the UVUDF at a 
depth comparable to the shallower ERS data (at which 
the UVUDF is highly complete), and expected number 
counts in the UVUDF to its full depth. We note that 
the O10 and H10 LFs included substantial corrections 
for incompleteness, so we do not expect the number of 
galaxies predicted in the ERS to match the observations. 

Predicted source counts for each selection method and 
dropout filter are presented in Table [3l In comparing 
the predictions to the observations, it is important to 
consider cosmic variance. For LBGs at these redshifts in 
fields the size of ERS and UDF, cosm ic variance could 
be a large effect (~20-30%, bias=1.5, i Somerville et al.l 
120041: IRafelski et al.ll2009HMoster et aLlboill ). However, 
in practice these fields are so close together on the sky 
(see Figure [3]) that they are not independent in terms of 
large scale structure. In addition, we do not correct for 
incompleteness effects at the faint end of the UVUDF 
data. With these caveats in mind, we conclude that the 
predictions are generally consistent with the LBGs that 
we observe. 

6.2. Resolving Galaxy Structure 

The deep WFC3 UVIS data provide the depth and 
resolution that allow us to study star-forming regions 
at z ~ 1 in the rest-frame UV. We identified a sam- 
ple of 179 galaxies with m 2 75 < 27.5 and 0.5 < 
Zphnt < 1-5, where p hotometric redshifts are taken from 
IRafelski et al.l (|2009t) . We find that galaxies frequently 
exhibit UV irregular morphologies and compact sizes 
(Figure H5j). with a median effective radius of 0'.'19±0'.'01 
(1.5 kpc at z = 1) in the F275W filter. The F275W sizes 
are broadly consistent with those measured at rest-frame 
~ 4000A which is probed by ACS I-band for 0.5 < z < 1 
and z-band for 1 < z < 1.5. At these wavelengths, we 
find a median size of 0.18 ± 0.01, suggesting that the 
recent star formation is occurring on the same spatial 
scale as previous generations of stars. However, when we 
deconstruct the galaxies clump-by-clump, clear morpho- 
logical differences begin to emerge. 

In Figure [16] we show some examples of z ~ 1 galax- 
ies. The object in the second row is at z p hot = 0.67, 
where the F275W probes the light from short-lived O 
and B stars. In the UV, most of the light is concen- 
trated in a few bright clumps. However, images at longer 
wavelengths (left column of Figure [TBI reveal that these 
clumps are within the disk of a well defined spiral galaxy, 




r„ (") 



Fig. 15. — Distribution of effective radii for galaxies in the 
UVUDF with m 2 75 < 27.5 and 0.5 < z phot < 1.5. We plot both 
the F275W (solid) and I-band (long-dashed) distributions, with 
the resolution limit marked with a vertical dashed line. Despite 
visible differences in the morphology between the rest-frame UV 
and rest-frame optical, the sizes remain approximately constant. 

with a clear bulge component at the center. If seen only 
in the UV, this object may res emble the clump-cluste r 
galaxies observed at z > 2 (e.g. Elmegreen et alj feOOS). 
and p redicted to form by fr agmentation within gas rich 
disks (jCeverino et al.l [2012). Clumps are predicted to 
migrate toward the center of the disk and coaelesce to 
form a bulge, whi ch eventually should stabilize the disk 
(jDekel et al.ll2009D . 

UV-bright clumps were previously seen in the same 
object using WFPC2/F300W images (Voyer et al 2009), 
but the significantly higher resolution of the WFC3 UVIS 
data (WFPC2/F300W FWHM=0'/27 compared to the 
WFC3/F275W FWHM=0'.'ll), enables us to measure 
star-forming regions as small as ~0.8 kpc (at z = 0.67), 
reaching 8a above the background level. One of the 
clumps that is unresolved in the WFPC2 images is clearly 
resolved into two clumps with diameters of 1.0 kpc and 
1.5 kpc in the WFC3 image. We have also identified 
clumps that do not appear to reside in a larger op- 
tical disk (bottom row Figure [T6|) . This object is at 
•Zphot = 1.18 and contains clumps with sizes ranging from 
0.7 -1.6 kpc. 

7. SUMMARY 

The UVUDF project obtained WFC3/UVIS observa- 
tions of the Hubble Ultradeep Field in three NUV fil- 
ters, F225W, F275W, and F336W (Figure [TJ). The UDF 
was observed with each filter for a total of thirty or- 
bits. The data were taken in three observing epochs 
with three orientation angles (Figure [2J . The data in 
the first two epochs were taken with 2x2 binning of 
the CCD readout in order to reduce the read noise that 
limited the sensitivity. For Epoch 3, as described in this 
paper (Section 13. 1|) . the observing strategy was changed 
to use the WFC3 post-flash capability to add additional 
background to the observations in order to mitigate the 
degredation of the CCD charge transfer efficiency. The 
post-flashed data were taken without binning the CCD 
readout, because the additional background noise domi- 
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Fig. 16. — HST gallery of clumpy spiral and irregular galaxies 
in the HUDF. (From left to right) Each panel shows (a) ACS BViz 
color combined image, UV WFPC2 F300W images from Voyer et 
al. (2009), and UV WFC3 F275W image of the combined Epoch 1 
and 2. All images are 5" x5". Photometric redshifts of each galaxy 
from top to bottom are 0.63, 0.57, 0.77, 1.18. 

nated the read noise. Coordinated parallel observations 
were obtained with ACS /WFC3 in order to provide very 
deep B-band fields, and the Epoch 3 parallels fall on top 
of one of the HUDF09 deep optical parallel fields (Fig- 
ure ©. 

The UVUDF observations present several data reduc- 
tion challenges. The team has produced new calibration 
files for the binned data using binned and unbinned cali- 
bration observations obtained by STScI. In addition, we 
have reprocessed darks for all the data, modeling each 
dark's gradient and using a more aggressive approach to 
flag hot and warm pixels (see Section 14. 1[) . 

The UVUDF data provide a stark demonstration of 
the effects of charge transfer inefficiency In this paper, 
we provide evidence that CTI causes increased scatter 
in the photometry of sources far from the CCD readout 
(Figures |H[5J and [SI- The application of a statistical cor- 
rection to the source flux densities based upon distance 
from readout is shown to reduce the scatter but not down 
to the level predicted from sky and read noise. We also 
find evidence that some faint sources far from the readout 
are lost to traps on the CCD in the data that were taken 
without post-flash (Figures [51 and [TUJ) . We agree with 
the STScI recommendation that future WFC3 UVIS ob- 
servations that require very sensitive measurements use 
the post-flash. The CTI also has demonstrated effects on 



the observed shapes of sources in the UVIS images, elon- 
gating them in the direction of the readout (Figure [7]) . 
This effect is problematic for both astrometric alignment 
and morphological analysis. STScI has released a prelim- 
inary version of software to apply a pixel-based correc- 
tion for the CTI, but it will need significant testing and 
verification before it is stable enough to justify its use in 
producing enhanced science products for the archive. 

The UVUDF data complete HST's panchromatic cov- 
erage of the Hubble Ultradeep Field. These data are 
applicable to a wide range of science topics. The mea- 
surement of the UV luminosity function, together with 
the mass function measured at longer wavelengths, will 
provide a statistical picture of the history of star forma- 
tion during its peak epoch. The superb spatial resolu- 
tion of UVIS will allow detailed analysis of star-forming 
"clumps" in galaxies, extending results obtained from op- 
tical images of z ~ 2 sources to later times and exploring 
the build up of normal galaxies. The UV sensitivity and 
spatial resolution will provide vital tests of the escape 
fraction of Lyman continuum photons from sources at 
z < 3, and of the star formation rate efficiency of neutral 
atomic-dominated hydrogen gas at z ~ 1 — 3. Finally, the 
new UV measurements enable significant improvements 
in the estimation of photometric redshifts. These several 
science investigations will be presented by the UVUDF 
team in later papers. 

In the current paper, we have presented preliminary 
results of the kinds that we can obtain. We used the 
UVUDF data to select Lyman break galaxies at red- 
shifts 1.7, 2.1, and 2.7. We find that the number density 
of dropouts is largely consistent with the number pre- 
dicted by the published luminosity functions based on 
measurements in the ERS. In addition, we confirm that 
UVUDF images of clumpy galaxies at z ~ 1 have suf- 
ficient sensitivity and spatial resolution to support the 
planned analysis of the evolution of star-forming clumps. 

There are many science uses for UV images of the UDF 
beyond those outlined above. This Treasury project will 
support a broad range of archival research. At the mo- 
ment, we are limited by the need to continue character- 
izing and correcting the CTI effects. We expect that the 
CTI correction software will become stable in the com- 
ing year. When it has been robustly verified, we will 
produce enhanced science products to be distributed by 
the Mikulski Archive for Space Telescopes (MAST). 

The UVUDF observations are currently the most sen- 
sitive component of a new generation of surveys that ex- 
ploit the unique capabilities of WFC3/UVIS. The sur- 
veys will leave a legacy of UV imaging for use in a wide 
range of research. The next logical step in expanding 
HST's UV legacy will be deep observations over a wider 
area than the UDF, in order to sample the variety of 
galaxy populations and their environments. These vital 
observations will greatly augment studies with the next 
generation of telescopes such as ALMA and JWST. 
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